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Abstract 

The relation between standard Loop Quantum Cosmology and full Loop Quantum Gravity fails already 
at the first nontrivial step: The configuration space of Loop Quantum Cosmology can not be embedded into 
the configuration space of full Loop Quantum Gravity due to a topological obstruction. We investigate this 
obstruction in detail, because many topological obstructions are the source of physical effects. For this we 
derive the topology of a large class of subspaces of the Loop Quantum Gravity configuration space. This 
allows us to find the extension of the standard Loop Quantum Cosmology configuration space that admits an 
embedding in agreement with [1] . We then construct the embedding for flat FRW Loop Quantum Cosmology 
and find that it coincides asymptotically with standard LQC. 
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1 Introduction 

The construction of a UV-complete theory of quantum gravity is one of the main open problems in theoretical 
physics, that sparked research programs in a variety of directions. A particularly fruitful program is Loop 
Quantum Gravity (LQG) [2, 3, 4], as it represents a mathematically well developed quantum field theoretic 
framework. An important source of insight for the development of LQG is Loop Quantum Cosmology (LQC) 
which is a symmetry reduced quantum mechanical model that exhibits many features of full LQG. The expec- 
tation that LQC captures essential features of LQG is often motivated by the following argument: On the one 
hand one expects that quantum effects are most important near a singularity of classical General Relativity 
(GR); on the other hand one expects in light of the BKL conjecture [5, 6, 7] that the dynamics near a singu- 
larity of GR is well approximated by the dynamics of decoupled homogeneous cosmologies, so one may expect 
that LQC provides essential insight for the dynamics of full LQG in situations where the quantum effects are 
expected to become most significant. 
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Symmetry reduced models of LQG, in particular LQC, are obtained as "loop quantizations" of symmetry 
reductions of the classical theory underlying LQG. This "classical reduction then quantization" procedure how- 
ever weakens the link between the reduced and full theory, because important features of the full quantum field 
theory may be overlooked by going through a classical symmetry reduction and re-quantization. Ideally, one 

A phenomenological investigation of this issue has been avoided in the LQC framework and the Spin Foam 
Cosmology (SFC) framework (see e.g. [10, 11]). It can be easily understood why SFC are insensitive: standard 
Spin Foams are not a path integral version of LQG. They are from the point of view of canonical LQG a 
restriction of full LQG to piece-wise linear LQG with low-valent vertices and unknotted edges. However, 
standard LQC is cmbeddable into full piecewise linear LQG 1 . The embeddability issue is thus avoided at the 
very start of the program. This is the technical reason why the investigation in this paper is independent of the 
significant advances in SPC. 

Canonical LQC is on the other hand sensitive to the issue. It seems that the standard LQC avoided the 
issue so far because only the topological obstruction to embeddability is known and the work-around in terms of 
piecewise linear LQG exists (we will comment on this at the beginning of section 5.2); however exploring possible 
deviations from standard LQC should be worthwhile. It is the purpose of this paper perform the first step into 
this direction. For this it is first necessary to investigate the topological origin of the obstruction to embedding 
standard LQC into LQG in detail and then to provide an explicit construction of an extension of standard 
LQC that permits an embedding. The investigation of this mathematical problem is physically motivated by 
the fact that configuration space topology is the source of interesting effects e.g. in solid state physics or in 
Euclidean QFT. Moreover, one of the most celebrated results of full LQG is the kinematic discreteness of the 
area operator. This can be understood as a direct consequence of the topology of the configuration space. We 
can in light of these precedences not a priori exclude that an embeddable version of LQC admits a very different 
phenomenology than standard LQC. The investigation of this phenomenology is however a vast subject and is 
thus beyond the scope of this paper. The main results of this paper are 

1. The induced topology on one-dimensional affine subspaces of the configuration space of full LQG is the 
spectrum of the algebra of continuous asymptotically almost periodic functions. This topology is finer than 
the Bohr compactification of the group (R, +), which is used to model the compactness of the configuration 
space of full LQG in standard LQC. 

2. We construct an embeddable version of flat FRW LQC, which differs from standard LQC by the intro- 
duction of configuration operators that vanish at infinity. However, standard LQC and embeddable LQC 
coincide asymptotically. Standard LQC thus captures the universal holonomy modifications. 

3. We do not see any obstruction to generalizing our construction of an embeddable flat FRW LQC to other 
models such as [12]. 

The paper is structured as follows: We first perform necessary asymptotic analysis in sections 2-4, which can 
be skipped by a physically interested reader. The content of the individual sections is as follows: 

In section 2 we reexamine the issue of non-embcddability of the standard LQC configurations space into 
the configuration space of full LQG. As found in [8], the obstruction to embedability of standard (Bianchi 
I, isotropic) LQC is the violation of almost periodic dependence of general LQG spinnetwork functions when 
evaluated on the minisupcrspace variables used in LQC. As it turns out [1] the solution is to extend the standard 
configuration space of LQC [13] to cylindrical functions, which are supported on arbitrary edges, not only on 
edges which are straight with respect to the metric of a chosen cosmological background. From the mathematical 
side the problem of constructing this extension is equivalent to solving the differential equation for the parallel 
transport (in this context also referred to as holonomy) of the Ashtekar connection along arbitrary edges. As 
usual on minisupcrspace the Ashtekar connection can be parametrized by a real parameter c. As a warm up we 
apply the first obvious approach. That is, we apply slight perturbations to straight edges. This leads to a series 
expansion of the general solution into powers of a perturbation parameter e, where by construction the zeroth 
order is standard LQC construction. Demanding the perturbation contributions to remain small for arbitrary 
large values of c turns out to require e ~ c _1 . Hence the perturbation approach already hints to a solution 
where contributions of perturbations lead to corrections as inverse powers of c. 

1 Piecewise linear LQC [9] is an interesting modification of standard Loop Quantum Gravity, where arbitrary graphs are replaced 
by piecewise linear graphs, where edges can be knotted and vertices can possess arbitrarily high valence. The kinematic theory 
differs significantly form standard Loop Quantum Gravity, while the diffcomorphism-invariant theory can be shown to be equivalent 
to diffcomorphism-invariant standard LQG if a particular extension of the diffcomorphism group is chosen. However, the precise 
construction of the path groupoid and extension of the diffcomorphism group may have observable consequences and it is thus not 
for a theorist but for experiment to decide which construction if any is realized in nature. 
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In section 3 we use the Liouvillc-Grecn-Ansatz [14] to determine the the dependence of an arbitrary spin 
network function on the minisuperspace variables. For this the general solution to the holonomy ODE in terms 
of a series expansion in inverse powers of c is constructed to arbitrary finite order in c _1 . We provide an explicit 
finite upper bound for the error of truncating the series at arbitrary finite order. 

In section 4 the limit c — > oo of the constructed solution is analyzed. It is shown to coincide with the 
standard LQC construction. We show that the desired extension of the LQC configuration space consists of the 
standard part of functions almost periodic in the parameter c plus functions which vanish for c = and c — > oo 
as observed in [1]. 

In section 5 we use this result, to construct an explicit embedding of isotropic Bianchi I LQC into LQG and 
discuss its relation with standard LQC. 

In section 6 we conclude with remarks on the physical interpretation of this result. We close our presentation 
in section 6.1 with an outlook on future work. 

To complete our presentation, the appendix gives parts of the explicit computations. Additionally an 
alternative derivation for the geometric interpretation of parameter c in terms of scalar curvature is presented, 
which uses a recently developed coordinate- free description [15] for the Ashtekar variables. 



2 Setup 

In this section we give a brief introduction to the symmetric setup used for the construction of standard LQC 
[13]. In appendix C a more detailed coordinate free treatment due to [15] is provided. 



2.1 Holonomy ODE for Homogeneous Isotropic Cosmological Model 

Assume a 3 + 1-foliation of a globally hyperbolic four dimensional space time Ai <~ R x E, where S denotes 3- 
dimensional spatial Cauchy surfaces. Let an edge e C E be given. Assume a trivial principal SU(2) fibre bundle 
and choose a coordinate chart covering a subset of E containing e. Then then pull-back of the real Ashtekar- 
Barbero connection as a su(2)-valued one form to E can locally be written as A(x) = A^x) dx a (g> tj, where 
{t/}/ = i j 2,3 denotes a basis of su(2). The embedded edge e can be written as a map e : R D [0, T] 3 1 1-> e(t) C E. 
The parallel transport of A along the edge e is referred to as its holonomy and is defined by the following ODE 
(we follow the conventions of [8]). 

^_h(e(t)) = -A(e{t)) h(e(t)) with initial condition h(e(0)) = l S u(2) (2.1) 

and A(e(i)) = A ! a {e{t)) e a (t) ® tj. In what follows we will denote derivatives with respect to t by dots (e.g. 
e(i) := (e a d a )(t) = g|e(t). Also we frequently suppress the dependence on t to shorten our notation. Hence we 
write for (2.1) 

h = -A(e)h . 

In this paper, we will be particularly concerned with homogeneous cosmology. This is obtained through a simple 
transitive action of a 3-dimcnsional Lie- group G on E, which allows the identification of E with G. Using this 
identification an invariant basis {x } a= i j2 ,3 in TE can be described by left / right-invariant vector fields on G. 
If one follows an integral curve X of one of these vector fields on G respectively E, then the components of the 
metric tensor are constant along that curve: 

(x a ,x b )\ x = g ab \ x = const 

Given the invariant basis {x a } a=1 , 2 ,3, the dual invariant basis of T*E, {x b } b=1 . 2 ,3, can be obtained from its 
definition x b (x a ) = 5 h a . For Bianchi I, we use G = M 3 . We can furthermore impose isotropy by enlarging the 
isometry group to E(3). Using a general result on symmetric connections [16], one obtains a parametrization 
of isotropic Bianchi I connections by A = c ■ 5 ! a dx a ® rj. Using e = e a (t) d a we find A(e) — c • 5^ e a ® r/. Using 
the defining representation of 51/(2), the holonomy h can be written in matrix form 



—b a 



with |a| 2 + |6| 2 = 1 
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Taking the usual basis of su(2) 

I 

Ti = i 



V 





■1 





( o ^ 






°) 


t 2 = i 




r 3 = i 






/ 






v ° 


V 



and using the shorthands 



e 1 ^) =: x 
e 2 (t) =: y 
e 3 (t) =: z 



m — x — ly 



77 = Z 



we can write (2.1) as: 



a b 



V 



-b a 



n m 



-ic 



m — n 



V 



a b 
-b a 



with initial conditions: a(0) = 1 , 6(0) = . From that we obtain two first order ODE's 

a = i c (na — mb) 
b = i c (nb + ma) . 

This can be transformed into a second order ODE for a. 

a = N a + M a with M := ^ 

m 

and TV := ic(n — Mn + ic(mm + n 2 )) 

where the initial condition for a now read as: 

a(0) = 1 d(0) = ic n(0) 



Using the transformation 
we rewrite (2.4) to 



d = 



d 



m/0 



[\m 2 - ^m + n) d 



(2.2) 



(2.3) 
(2.4) 

(2.5) 
(2.6) 

(2.7) 



.4 2 

Equation (2.7) is a linear ODE of second order. Its solution d(t) consists of a linear combination of two 
fundamental solutions d^(t) and d^~\t), that is 

d = d< + >+d<-> . 

With this, the solution to (2.4) can be written as 

o= Vm~ 04 ( +)d (+) > ( 2 - 8 ) 

where ^4(+), ^4(_) are constants to be chosen such that the initial conditions (2.5) are satisfied. 

2.2 Non-Embedability of Configuration Spaces 

In [8] it was shown, that the configurations space of LQC, given by RboIh-j the Bohr compactification of the 
real line cannot be continuously embedded into the configuration space of LQG, the space A of generalized 
connections. It was found that in order to retain continuity it is necessary that the solutions (2.8) depend 
almost periodic on the parameter c for arbitrary edges, not only edges which are straight with respect to the 
symmetric background. However it was shown explicitly that generically this is not the case. 
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Spiral Arcs. In [8] an exact solution to the ODE (2.5) it was found for spiral edges. This is the most 
general case in which the ODE has constant coefficients, because in the spiral case M(i) = M(0) =: Mo and 
n(f) = n(0) =: uq are constants. Anticipating the notation of section 3 the ODE (3.3) reads 



d = n 2 A 2 d with A 2 =A 2 ( K ):=l + ^ + % 



(2.9) 



where k :— i c, cto := —M Q no and po := jMq and we have used arc-length parametrization |m| 2 + n 2 = 1. It has 
two fundamental solutions d^{n, t) — <£ iKAt and we can obtain a solution to the holonomy ODE (2.5) and to the 
original ODE system (2.3) by analogy with section 3.2 by imposing initial conditions (a(n, 0) = 1, cl(k, 0) = Kno, 
b(n, 0) = 0, b(n, 0) = klto(O)). The solution 2 is given by 



/ % 1 m(t) 



E 

?=±i 



i + e 



2 k no — Mo 
2kA 



b(M) = y/m(t)m(0) J2 { l -^ KtK } 
e=±i ^ J 



(2.10) 



To obtain an explicit parametrization for the solutions (2.10), we set m(t) = /z-<e (w.l.o.g. u € (D, A e R), hence 



Mo = iA and m(0) = /i. Moreover let no := v with veE, |//| ; 



1 and A := Vc 2 A 2 = + c(c - i^A)] 



1/2 




Figure 1: Spiral arc as described by our setup. Recall that m = x — iy and n — z. We have R 2 = x 2 + y 2 = 
and the step height per period is given by D — 2 ^. The limit A — » corresponds to a line in the (a;,y)-plane 
with direction (3?(/i), — 5(/i)). The limit A — > oo corresponds to a line in z-direction. In case of a planar circle 
(y = 0) we have R 2 = A~ 2 . 

Then (2.10) reads 



a(K,t) = e 2 
b{K,t) 



v J cos(At) + ^(2^ c - A) sin (At) 



7 A/ 



1C 



/L* € 2 "" ^ sin(At) 



(2.11) 



Obviously (2.11) exhibits a non-almost dependence of o(k, t), b(n,t) on c, unless the underlying curve is a line 
(in this case A = and v = 0). However one observes from (2.11) that in the limit c — > oo periodicity is 
asymptotically restored. This observation is the starting point for looking at the properties of the holonomies 
along arbitrary edges. The first obvious Ansatz is to look at the effect of small perturbations of straight edges. 



Perturbation of Straight Edges The detailed computation can be found in section B. Here we only quote 
the result. The general setup is given as follows: Let eo(i) = (e 1 (i), e 2 (t), eo 3 (t)) =: (xq, yo, z o) ■ Now 
assume the edge 70 is deformed into another edge 7 such that 7(t) = (xq + ex, yo + ey, zo + ez), where 
x = x(t), y — y(t), z = z(t), and e — const is a small deformation parameter. For eo being a straight line 
and under some simplifying assumptions described in section B.5 one obtains a solution to (2.8) in terms of a 
formal power series in the perturbation parameter e. However one finds from the computation in section B that 

2 We denote it by German type letters. 



G 



in order to ensure "small" contributions from the perturbation even for large c one has to require that e ~ c 1 . 
One then obtains a solution to (2.8) which looks as follows (see section B.5 for computational details): 

o(/s,t) - J2 «E ,T,t *{l + 0(c- 1 )} • (2.12) 

<T = ± 

From (2.12) it becomes obvious that the property of asymptotic almost periodicity carried by the solutions 
(2.10) is not a special property of considering spiral arcs. Rather it is a generic property valid for general 
perturbations of straight edges. Hence we are referred to look for a general solution of (2.10) for large c in terms 
of an inverse power series of c. Fortunately, there is well developed mathematical framework available (see e.g. 
[14]), which uses the so called "Liouville Green Ansatz". This is employed in the next section. 



3 Solution by a Liouville Green Ansatz 
3.1 Liouville Green Ansatz 

In this section we are going to compute an expression for the holonomy along arbitrary edges. It is based on 
the presentation in [14]. In principle, this approach is viable for general differential equations of the form 

a(K,t) = n 2n f(n,t) a(K,t) (3.1) 

where n G N, k e (D is a complex parameter and /(k, t) is an analytic function of t £ (D with uniform asymptotic 
expansion 

/(K,t) = / (t) + M) + ^W+... for |k|-xx>. (3.2) 
Using (2.4) we can rewrite equation (2.7) as 

d{n,t) = K 2 (\ + ^ + fWp\d{K,t) with 1 = f (t)=m{t)m(t)+n{t) 2 

a(t) := f 1 (t)=h(t)-M(t)n(t) 



m ■= Ht) = \M\t)-\M{t) 



(3.3) 



fk(t)=0 if 



k > 2 



and t € 1R, k = ic and we have chosen arc-length parametrization in order to set m(t) m(t) + n(t) 2 = 1 . For 
shortness of the notation we will from now on suppress the dependence on k, t where possible. 

Now the strategy is to write (3.3) as 

d-K 2 d=(na + P)d (3.4) 

and to regard the rhs as an inhomogeneity. Then two fundamental solutions d^ (n, t) to (3.4) can be obtained 
from the Liouville-Green- Ansatz 

d W {K ,t)=^X^\t)^^^j , (3.5) 

where £ = ±1 and A^\t) = const = 1. Plugging this Ansatz into (3.4) we obtain 
dW = <e^*{ K 2 AW + k(X^A { ( ] + 2£XW) + (X^A^ + 2£{X^A ( ( ) + X^A^) + JfW) 

oo 

+ E K ~ r ( xWA % + MX W 4 e li + + 2£X^A% + 2X^A^ + X^A™) } 

r=l 

and 

oo 

(« 2 + na + P) dW = e iKt X^{k 2 + k(A^ + a) + + A { 2 e) + ]T K ~ r (a% + aA% + PA™) } 
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(3.6) 



Now we compare coefficients for every order of k. 

0(k 2 ) : = 

0{n r ) : 2£X^ = aX® 

O(k ) : IW +2£A[ e) xW = ^jW 

These equations can be integrated. From ©(k 1 ) we find 3 

X( e \t) = X { e) «2 So <*{>)*> with X^aa integration constant . (3.7) 
yielding the following identities 

XW £ l(<) I l 

__ = _ a and __ = _a+-a . (3.8) 
We can use (3.8) in order to rewrite the O(kP) and C(K _r )-expressions of (3.6) to get 

4 £) = Uf>-U-±cA (3-9) 



r+1 2 \ 2 4. 
which can be integrated for r > to 



2 V 2 4 

4<> =£(/?-£ d - £ a 2 ) AW - - £ AW (3.10) 



4')(t) = cf ) + ££(/3( S )-£d( S )-i a 2 ( S )) (3.11) 

4%(t) = ci%-\a{t)AVXt)- l -AW{t) + l -j^^ (3.12) 

where C^C^ are integration constants. Here we have used the fact that £ 2 = 1. Notice the different signs 
in front of the a term in the integrals on the right hand side of (3.11), (3.12), which result from a partial 
integration of the aAi^ '-term in (3.10). 

The integration constants (7^,(7^ and X^ can be fixed by imposing the initial conditions (2.5) and 
using the Ansatz (2.8). 

3.2 Power-Series Solution to Holonomy ODE for General Edges 

We can now construct a formal solution to the original ODE (2.5) and to the original ODE system (2.3). To 
summarize, we have 

a(n,t) = y/m{t) (A (+) d (+) (K,i)+A ( _ ) d ( - ) (K,i)) (3.13) 
a(n,t) = £ M (t) o(k, t) + y/m(t) (A {+) d {+) (K,t) + A { _ } S~\K,t)) . (3.14) 



with initial conditions 
and respectively 

with initial conditions 

and (k, t) given by (3.5) 



a(/c,0) = l d(/c,0) = /cn(0) (3.15) 

b(n,t) = yjm(t) (B (+) d (+ \n,t) +B { _ ) (f- ) (K,t)) (3.16) 
6(k,0)=0 &(/c,0) = /cm(0) (3.17) 



3 Here and in what follows we will always integrate from to t and insert the appropriate integration constants when necessary. 
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Fixing Xq\ Without loss of generality we can set the integration constant X^p in (3.7) to — 1, because 
as a constant it can be absorbed into A-(-) respectively £(+), -B(-) in (3.13) and (3.16). Therefore we have 



3.2.1 Fixing integration constants and solution for b(n,t) 
We have 

0) = = v^O) £ ^ + B ( _, £ ^ 



- ^™ +) (o) , o ^4" } (o) 



n=0 



n=0 



- — S 



(-) 



Now we compare coefficients for every order of k. 

O{k ) : A ( e) = 1, hence = + B { _ } 

O(K- n )\ n>0 : With B (+) = -B { _) we have = A { + } (0) - ^(0) ~» Al +) (0) = Ai _) (0) 
Additionally we have 

6(/s,0) = /sm(0) - ^^6(k,0) + 7^(0) (s (+) ci (+) (K,0) + B ( _ ) d ( " ) (K,0)) 
(3 ^ 0) v^(0)B (+) (d< + )(/c,0)-d<-)(/t,0)) 



Comparison of coefficients gives 

Oin 1 ) : A ( ( p = 1, hence 

O(k ) : 

O(«- n )|n>0 : A { + ] =aV 
If we compare this to (3.11), (3.10) we get 

ci +) = c^ = -^ and c w =Ii(-) (0) , c^ 1 = -ii(+) ( o) 

3.2.2 Solution for a(/c,t) 
We have 



2 v /m(0)B (+) ~» 
a(0) + 2 4 +) (0) ~» 


B {+) ~ B (-) = 2 




4+)( 0) = 4- ) (o) = -^ 


4"li(o) = - ^ ^(0) - ii^ +) (o) + M^W 



«(«, 0) = 1 = v^(0) L +) £ ^1 + A { _ } £ 



° Ai~\0) 



n=0 



n=0 



Now we compare coefficients for every order of n. 

O(k ) : A^=l, hence 1 = v / ™(0)(^(+) + A (-)) 



O(«-")|n>0 = 



= A (+) 4 +) (0) + V) A ™ _) ( ) ~> 



= [m(0)]"3 - 



i-»(o) = -^i + »(o) 



(3.19) 



(3.20) 



L n=0 n=0 n=0 J 



(3.22) 



(3.23) 



(3.24) 



(3.25) 
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Additionally we have 

d(K, 0) = Kffl(O) 



M(0) 



a(n,0) + v^(0) cZ (+ >(k,0) + d (_) (K,0)) 



(3.24) 
(3.25) 



M(0) 



( (0)\y j +) (0) | y^ il +) (0) 



n=0 



n=0 



-4 



(-) 



Es±n yo) \ ^ /In (> 
AC™ ^ K™ 

n=0 n=0 



(3.26) 



Comparison of coefficients gives (A^ = 1) 



©(k 1 ) : n(0) = v^0)(A (+) -A ( _)) 



(3.25) 



_ l±rt(0) 



2 A /m(0) 



O(k ) 



0(«T n ) 



= 



M(0) 
2 



+ v^(o){V)(4 +) (o) + 2M) - A-)(4 _) (o) + ^)} 



^WffVi - a(0)"(0)+M(0) 
A l W _ 2(l±n(0)) 



|n>0 



= A, 



(+) 



4+i(o) + £ ^ I 4 +) (o) + 4 +) (o) 



Ai-^(0) + ^Ai-\0)-A n -\0) 



(3.27) 



(3.25) 
^> 




-\(0) = 


"(0) 4 M 


- } (o)- 


2 A « W 2i (+) /1 » 


>(0) 


(3.25) 




- } i(0) = 


a(0)^(- 


~ } (0) + 


ii(-)(0) + i^il + 


} (0) 



where we have just multiplied the second to last row by —^47-7 in order to arrive at the last row. If we compare 
this to (3.11), (3.10) for t = we get 



(-) 



(+) 



*(0) n(0) + M(0) 
2(n(0) + l) 



and cff 1 = -i4 tl 4- ) (0)= ^^^^(0) 
n+1 2 " w 2 n(0) + 1 



( -)_ a (0)»(0) + M(Q) ( _)_ U( + ) i(+lf] _ l n(0) + l 1(+) 

Cl _ 2(n(0)-l) 2i H in (U) " 2n(0)-l " (U) 



(-) 

3.2.3 Final Solution 

We have thus constructed a formal solution to the ODE-system (2.3) with initial conditions (2.5) 



a(n,t) = y/m(t) (A (+) d i+) (n,t) + A^ d { ~ } (n,t)) 
b(K,t) = (B {+) d(+\K,t) + B { _ } S-\K,t)) . 

Using arc length parametrization m(t) m(t) + n(t) 2 = 1 and 



a(t) := h(t) - M(t) n(t) (3(t) := ^M 2 (t) - ^M(t) 



we then have for I = ±1 



\r=0 



with 



(3.29) 



(3.30) 



A%) = <%> + {fl(^{s)-^{s)-\a\s)Y^{s)ds 

4%® = cW 1 -\a{t)A¥Xt)- i -AW{t)+ t -J^^ , (3.31) 
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where the last line is valid for r > 1. For a(n,t) we have the integration constants 

n(0)+l (/) _ a(0)n(0) + M(0) w _ £ n(0) - £ - ( _ <} 

w V^gy 1 " 2(n(o)+£) 6r+l _ 2 ' y^yr^ (0) 

and for 6(k, t) we have the integration constants 

B W =l-^ cf> = -« ^ = 1.^(0) . 

Given these solutions in terms of a formal power series it remains to discuss the finiteness properties of this 
series. This issue will be discussed in the next section. 



4 Asymptotics and Finiteness of the Solution 

4.1 Asymptotic Behavior of SU(2) Holonomies 

Using the Liouville- Green method and Horn asymptotics (see e.g. [14]) it is straight forward to find the 
asymptotic behavior of the holonomy differential equation. We present the general procedure here and refer to 
appendix A for details. 

Recall that the general second order linear ordinary differential equation a(t) = A(t)a{t) + B(t)a(t) is 
transformed into standard form d = (B(t) + (jA 2 (t) - ±A(t)))d(t) using a(t) = exp(| f* dsA(s)). This 
transforms the holonomy differential equation (in arc-length parametrization for t) into 

d{t)={n 2 + na{t)+p{t))d{t), (4.1) 

where we introduced n := ic. The Liouvillc-Grccn Ansatz for the solution is 



d^(t) = exp ±nt±U ds a{s)\ 1 + V ^-J!Z , (4.2) 




which suggests that d^(t) — > C^e Kt as |k| — > oo, with prcfactor To prove this behavior, one needs 

convergence of this Ansatz, which is rarely given. In fact the Liouville-Green Ansatz for a generic holonomy 

diverges. However, any truncated Ansatz X)fc=i Ak K ^ i s wen behaved (if the n-th derivatives of a, (3 are 
bounded) and hence we consider the truncation error 

4 ±) (*) = ^ (±) w-E^#- ( 4 - 3 ) 

fc=l K 

The asymptotic behavior can the be obtained from bounding the truncation error using the following steps 

1. Insert the truncation error into the ODE to obtain an equation for e^(i). This ODE can be solved using 
variation of constants in Horn's asymptotic Ansatz Z ± = (k+ \<x) 2 exp (±(k£ + \ J a)). 

2. The key of Horn's Ansatz is that it solves the ODE to all positive powers of k, hence the kernel of the 
integral equation obtained form the variation of constants does not contain any positive powers of k and 
thus yields finite bounds on the solution when \n\ — > oo. 

3. Finiteness of the coefficients of the truncated Liouville-Green series together with the bound on the 
truncation error yields the anticipated asymptotic behavior of d^(t) C^e ±Kt as \k\ — > oo and thus 
also for a(t) = d(t) exp(i J* dsA(s)), since A is independent of k. 

4. The asymptotic behavior of the matrix elements of the holonomy is then determined by the asymptotic 
behavior of the linear combination of d^\t) that satisfies the initial condition for the specific matrix 
element, which yields the anticipated asymptotic periodicity of the holonomy. 
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A significant simplification to estimate of the truncation error can be obtained by using methods developed in 
[17]. There it is assumed that the edges along which the holonomy is integrated are holomorphic 4 in a finite 
radius in the complex plane around the part of the real i-axes that is integrated, because this yields a very 
simple bound on the Liouville- Green coefficients {t). This method is applicable if one is able to satisfy initial 
conditions s.t. all integration constants in the Liouville Green recursion (3.31) vanish. Unfortunately this 
requirement is in general not compatible with the initial conditions (2.5). A solution to (2.3) fulfilling the initial 
condition (2.5) docs in general not have vanishing ci^ . This is only possible in very special cases 5 . 

Let us now discuss this result: By construction the series (3.30) provides a formal solution to the ODE (3.4). 
However it is not not obvious whether this series is finite nor whether it converges. 



Finiteness. As explicitly shown in appendix section A, for analytic edges the series (3.30) is finite: at every 
finite inverse order of k the rest term 

e^ («,t) := dWlis.iJ-e^JJ-W^^^f (4.4) 



K 
k=0 



obeys the explicit bound 



, m , ex PU ms)\as / ft \ / (\ ft 

\ei e) (n,t)\ < ^-j '-^ |AW( K ,r)|dTjexp(exp^- ^ \a(s)\ds 

2$(t) 2 J^(r) \A^\r)\dr / $(*) g | X («, s)\ ds ' 

\(n,t) 6XP 1 \{ K ,t) 



X(K,t) 



where A(«,t) := inf < r < t {k+^(r)|} and \A^(k, t)\ = \e lKt X^ (t) 2 -^§^\ < H"" 1 exp (/„' \a(s)\ ds) \A^(t)\. 
Here we have set X^(t) = <s^fo a ( s ) ds according to section to section 3. The quantity x( K i s ) m given 
by x(k,s) = /3(s) - - 16(K 3 "i ( a S ( s))2 + 4 ( K +i S i( s) ) and thc functions a(i),0(t) are given as in (3.3). In 

the last line we have introduced the shorthand <&(t) := exp (j^ J Q * |a(s)|ds^. Hence for fixed t we find that 

|e^(/c,t)| <~ 0{n~ n ). At the same time, given thc finiteness of the interval [0,i] C 1R and assuming the analy- 
ticity properties of the functions a(t),(3(t) as in section A. 2, it is obvious from (3.31) that for every r < oo it 
holds that jA^^)] is bounded from above, because it is a finite combination of lower order terms. 

Convergence. In section A.l a general estimate is given for ^^(t)] is derived which only depends on r. 
There the results of [17] are applied. These results only hold under quite restrictive assumptions on the recursion 
relation (3.31), in particular it requires the integration constants Cn^ to be identical zero. Moreover the functions 
a(t),j3(t) are required to be holomorphic in a small strip around the R-axis in thc complex plane, in order to 
use Cauchy's integral formula to replace derivatives with respect to t by contour integrals. Unfortunately these 
requirements are too restrictive to be applied for the ^^(t)) for general curves 7. Interestingly the requirement 
C [ a ] =0 Vn > 1 is fulfilled for spiral arcs, for which however the general explicit solution is known [8]. 



4.2 Asymptotic Behavior of Holonomy Matrix Elements for Symmetric Connec- 
tions 

In summary, the last subsection has shown: 

Lemma 1 Given the analyticity condition (A. 28), (A. 30) and quasi- arc-length parametrization 
n 2 (t) + \m(t)\ 2 — 1 there exists a t 2 s.t. the SU(2) holonomy matrix elements satisfy for t < t 2 : 



a(n,t) -> 4 / m ^?2 { cosh ( nt + j- f a(s) ds j - in(0) sinh ( nt + i f a(s)ds 
V m (°) I \ 2 Jo J V 2 Jo 

b(n,t) — > —iy/m(t) m(0) sinh (^nt + ^ J a(s) ds 



(4.5) 



4 In this case one can apply Cauchy's integral formula in order to replace derivatives of a(t), /3(t) in (3.31) by integral expressions. 
5 See the paragraph on Convergence at the end of this section. 
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in the limit |c| — > oo. 

Notice that this gives the asymptotics only. In order to rc-dcrive the property that |a(«, t\ 2 + \b(n, t)\ 2 = 1 one 
has to start from the general solutions of section 3.2.3, compute |a(/c, t)\ 2 and \b(n, i)| 2 and then take the limit 6 . 

An edge e is parametrized by < t < T, which allows us to define the inverse edge e _1 (t) := e(T — t). 
Lemma 1 thus implies that there exists t 2 (e) > and ^(e -1 ) > 0, where the asymptotics (4.5) holds. Wc can 
thus restrict our attention to the compact interval t 2 (e)/2 <t<T — t 2 {e~ 1 )/2. For any point t in this interval 
we are able to onsider e^~(A) := e(t + A) and e^~(A) := e(t — A), and lemma 1 states that there are t 2 , t 2 > 
such that the asymptotics (4.5) holds for < t < t \ on e/" and analogously for e^~, implying that any point 
t in the interval t 2 (e)/2 < t < T — t 2 (e _1 )/2 has an open neighborhood where the asymptotics holds. Hence 
by compactness of the interval t 2 (e)/2 < t < T — t 2 (e _1 )/2 one can find a finite open covering of this interval 
where the asymptotics holds for each element of the covering. Thus the asymptotics of the holonomy is given 
by a finite number of matrix products the matrix elements of each having the asymptotics (4.5), which imples 
that the holonomy of the entire edge has exponential asymptotics in k. 

Now we introduce: 

Definition 1 A complex-valued function is called asymptotically almost periodic, iff it can be written as the 
sum of a continuous almost periodic function and a continuous function that vanishes at infinity and at zero. 

Lemma 2 The continuous asymptotically almost periodic functions form an algebra over C, and each element 
splits uniquely into the almost periodic functions plus functions vanishing at infinity. 

proof: Let /j be asymptotically almost periodic, so there exist continuous almost periodic functions dj and 
continuous functions 6, vanishing at infinity, such that fa = ai + bi. 

algebra: af is asymptotically almost periodic, since for a e C: aa is almost periodic and ab vanishes at 
infinity, /i + /i is asymptotically almost periodic since a\ + a 2 is almost periodic and bi + 6 2 vanishes at infin- 
ity. /i/ 2 is asymptotically almost periodic, since aia 2 is almost periodic and ai& 2 +a 2 6i+6i& 2 vanishes at infinity. 

uniqueness: Assume / is asymptotically almost periodic and there exists continuous a\ ^ a 2 almost peri- 
odic and continuous b\ ^ & 2 s.t. / = cti + 6i = a 2 + 6 2 , so a\ — a 2 = 6 2 — b\, which implies that a\ — a 2 vanishes 
at infinity, hence a\ = a 2 , thus / i-)- a\ + &i is unique. If c = the original ODE system (2.3) has constant 
solutions which are almost periodic. □ 

Let us recall that a gauge-variant spin-network function T 7 is a finite collection of edges 7 = (ei, ...,e„) with 
a matrix-element of an irreducible representation p of the gauge group (SU (2) in the present case) associated 
with each edge, such that T^{A) = ]X =1 p 3i (h ei (A)) mini . Using that a matrix element of an irreducible repre- 
sentation is a polynomial function of an the matrix elements of the fundamental representation, we find using 
lemma 2: 

Theorem 4.1 Let cA* denote the symmetric connection as described in section 2.1. Given a spin network 
function T 7 s.t. each edge e, 6 7 can be written ei — e^.i o ... o a.^ such that each eij satisfies the analyticity 
condition (A. 28), then T 7 (c^4*) is an asymptotically almost periodic function of c. 

This brings up the idea, to group the set of all spin network functions T 7 into equivalence classes, where 
T 7 ~ T 1 , are called equivalent if T~ / (A)\ c a„ = T', (A)\ c a, , that is if they coincide when evaluated on a symmetric 
connection cA». In section 5.2 this observation will be used in order to construct an embedding of the extended 
configuration space of loop quantum cosmology (given by the set AAP(]& of asymptotically almost periodic 
functions) into Cyl(A), the configuration space of full loop quantum gravity. 

4.3 Generalization to the Full Configuration Space 

Remarkably, the previous construction generalizes straightforwardly to any 7 bounded reference connection A» : 
let us consider a one-dimensional subspace of A, the space of generalized connections. Let this subspace contain 

6 There are mixing terms which contribute then. 
7 Not necessarily symmetric. 
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connections of the form cA*, where A* is a bounded su(2)-valued reference one-form, i.e. its components (^4*)^ 
are bounded in a trivialization A* = {A^Y a dx a (g> tj . Using the notation 



n(t) = (A*f a (e(t))e"(t) 
m(t) = ((A*)l(e(t))-i(A*) 2 a (e(t))))e»(t), 

we impose the "quasi arc- length" parametrization condition n 2 (t) + |m(i)| 2 = 1, such that the holonomy ODE 
(2.1) becomes 

w(t) = M(t)w(t) + (ic(h(t) ~ M(t)n{t)) - c 2 ) w(t), (4.7) 

where M(t) = and w — a or b with initial condition a(0) = 1, b(0) — 0. We thus have a tool to investigate 
the compactification of the configuration space of LQG in a general "radial 8 "through the induced Gel'fand 
topology on the radial coordinate c. To re-phrase this: every gauge variant spin-network function T 7 exhibits 
an asymptotically almost periodic dependence on c if the connection A is parametrized 9 as cA» . 



5 Quantum Symmetry Reduction 
5.1 Symmetric Quantum Connections 

A classical connection to (on a trivial bundle) is determined by the connection one-form A = <t*uj obtained as 
its pullback under a section a. A connection one-form then defines a homeomorphism from the groupoid V of 
piecewise analytic paths to the gauge group G by 



A : V -> G , e i y h e (A), (5-1) 

where h e denotes the holonomy along the path e. A quantum connection is thus conveniently defined as an 
element of Hom^, G), which imposes no further restriction on the homeomorphism. This is in contrast to 
a classical connection where differentiability allows to reconstruct the connection components unambiguously 
form a much smaller subgroupoid V a , that is local and can be constructed using in each chart (U d , {<j) a } d =1 ) as 
follows: 

The base space is vi°^ = U. At each point x Q e U d we consider the d linearly independent coordinate 
directions d(j) a and the d corresponding families of integral curves e" t : {0,t) — > U d which we consider as arc- 
length parametrized in the Euclidean metric on R d . The elements of the groupoid V then consist of all finite 
concatenations of these integral curves modulus zero paths and the source map is s(e° t ) = x Q and range map 
is r(e° Co t ) = e% o t {t). The components of a classical connection are then by construction completely determined 
in terms of the restriction to V through the holonomy differential equation 

^(.t )t 7 = d t h el (A) . (5.2) 

Let us now consider the subspace S of classical connections that are invariant under the action of a symmetry 
group H. We can use equation 5.2 to define a classical symmetric connection (tautologically) as follows: 

Definition 2 (Classical Symmetry) A morphism from A : V — >• G is called classically M-symmetric iff there 
exists a diffeomorphism 4> and gauge transformation g and an element B e S s.t. A(e) = g^ 1 (s(e))h^ e - ) (B)g(r(e)) 
for all e € V . 

Since V is a subgroupoid of V, we can immediately apply this definition to quantum connections. However, 
using this definition does not seem to capture a useful notion of symmetry for quantum connections as can be 
seen from a simple example: let B S S and B ^ S and define the homeomorphism as the extension by groupoid 
composition of a generating set that contains a generating set of V Q and generators that can not be decomposed 
into elements of V a : 

h e (B) : eeV 

A:e^{ (5.3) 
h e (B) : eeV 



8 Wc understand radial in the sense of a fictitious Banach-space. 

9 For a general connection A £ A. with bounded components A ! a (x) < oo £ S and Va, / = 1,2,3 in a chosen trivialization, 
choose e.g. c := sup^ a I \A I a {x)\ and introduce (A t )^(x) := c~ 1 A I a (x). 



14 



where V Q denotes all elements of V which do not admit nontrivial decomposition that contains an element of 
V . This quantum connection does not appear symmetric for almost all paths, except for the ones in P a . We 
are thus lead to use a stronger definition of symmetry for quantum connections: 

Definition 3 (Quantum Symmetry) A morphism A : V — > G is called quantum H- symmetric iff there exists 
a diffeomorphism <ft and gauge transformation g and an element B G S s.t. A(e) = g~ 1 (s{e))h^ e }(B)g(r(e)) 
for all e eP. 

These two definitions of symmetry applied to the morphism defined by the holonomy of a classical connection 
coincide due to the invertability of 5.2, but weed out strange quantum connections as the one defined through 
equation 5.3. 

Notice that usual Loop Quantum Cosmology is motivated by a model of Loop Quantum Gravity defined on 
V , so the difference between these two definitions does not occur there and one usually uses definition 2. 

5.2 An Explicit Quantum Symmetry Reduction 

The configuration algebra of standard Loop Quantum Cosmology is the closed span of exponential functions 
of the connection parameter c, so the exponential functions are a dense (using the ||.| loo-norm) subset of the 
configuration operators. The assumed relation r between exponential functions e lct in standard Loop Quantum 
Cosmology with full Loop Quantum Gravity configuration operators is given by the observation that holonomy 
matrix ht elements along straight edges of length t exhibit this exponential dependence on c, so r(e lct ) = h t . 
This relation extends straightforwardly to the span of the exponential functions by linearity. To extend it to 
the closed span however, one needs to assume a bound ||K/)||lqg — ll/ll°o to ensure convergence. 

In this first attempt we set difficulties 10 concerning issues of convergence and completion aside for now, but 
caution that one might be forced to amend the definition of r to extend it to the closure. We will denote by 
C C (R) functions of compact support and by C (H) functions, which vanish at infinity. 

The construction strategy for the quantum embedding follows [18]: Let A(c) be the embedding of the 
symmetry reduced configuration space, parametrized by c into the space of all connections, so A can be used to 
pull-back any configuration observable (in particular cylindrical functions) T(A) to a function T(A(c)) on the 
symmetry reduced configuration space, which defines equivalence classes of configuration observables in the full 
theory whose pull-back to the symmetry reduced configuration space coincides. Thus, a symmetry reduction is 
encoded in these equivalence classes. To relate the symmetry reduced model with the full theory, one constructs 
an embedding of the symmetry reduced configuration observables /(c). For this we have to find a representative 
r(f) in the full theory for each equivalence class, such that K/)U( C ) = f( c )- 

The difference between our present treatment and standard loop quantum cosmology is that our treatment 
forces us not to restrict ourselves to continuous almost periodic functions f a G AP(M), but to continuous 
asymptotically almost periodic functions / G AAP(M). Given a continuous asymptotically almost periodic 
function / G AAP(R), we can (by definition) find a pair f a G AP(R), f a G C (R) s.t. f = f a + fo- Lemma 2 
provides that the pair f a , f Q is unique. In order to construct an explicit embedding of AAP(Wi) into Cyl(A), the 
configuration space of full LQG, we are now going to construct a generating set of functions, which allows us to 
describe every / G AAP(M) without referring to the explicit series (3.30). For this we start with the definition 
of some functions. 




_i i i i i i i i i i ^. 

A a—1 a b 6+1 B B+A B+a-l ■" X 



Figure 2: Functions (schematically) used in the construction of the embedding. 



"It is in particular not obvious whether the completion of the embedding in the full LQG-norm lies within the Hilbert space of 
full LQG. 
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For each a < b let us define Xa,b € C C (R) through 

: x < a- 1 
x — a + 1 : a — 1 < x < a 

1 : a < x < b 
b-x + 1 : b<x <b+l 

: x > b+ 1 

For each a < b and / e C C (IR) with supp(/) C [a, b] we introduce the continuous (b — a + 2)-periodic function 
fa.b € AP(R) defined through its values at the defining interval 



Xa,b( X ) :z 



(5.4) 



fa,b{x) := 



: a — 1 < x < a 
f(x) : a < x < b 
: b < x < b+ 1 



(5.5) 



and extended to all of R by (6 — a + 2)-periodicity 11 . For a specific positive function g(x) and each A < a— l,a < 
fo, 6 + 1 < B we define the continuous (£? — A) -periodic function G^' fe S e ^4P(M) defined on the initial period 

A < x < B as 



<f (*) ■■= < 



: A < x < a - 1 



9@) 



: a- 1 < x < 6+1 



(5.6) 



: fe+ 1 < x < B 



and extended to all of Rby (B — A)-periodicity. Let us furthermore assume g € C (R) s.t. Vx : g(x) > and 
3g > : g(x) < ffi for |a;| > x . 

Lemma 3 Every f € AAP(R) there is a Cauchy sequence of sums of exponential functions plus g times sums 
of exponential functions converging to f in ||.||oo- 

proof: We use that there is a unique split f = f a + fo with /„ e AP(M) and f a € C (R). Using the density of 
the exponential functions AP(R) we find a Cauchy sequence of finite sums of exponential functions converging 

to fa- 

For f we use that for every e > there exists a, b € M and / € C C (M) with supp(/) C [a, b] s.t. ||/ — 
/||oo < £• Using / = Xa,bfa,b, we find that for every e > there exists a, b and and f a ,b € AP(M) s.t. 
||/o— Xa.fc/a.bH < £• Moreover, using g(x) > and ^(x) < |f for x > x Q , we find for e > there exist A, £? e R s.t. 
II^G^if - Xa.bHoo < e- Hence from 

ll/o -SG^f/a.blloo = ||/o - (ffG^if -Xa,6 + Xo,6)/a,6|| < e (l + | |/a,fc | |oo) 

and 1 1 f a b | |oo < 1 1 fo | |oo + e and the density of the exponential functions in AP(M) we find that for every e' > 
there exists a sum of exponentials h yielding ||/ — gh\\oo < e'. ■ 

An Explicit Embedding. Now we use this result in order to construct an explicit embedding of AAP(R) 
into Cyl(A), the configuration space of full LQG. For this we use the spiral arc solution of section 2.2 and the 
defining representation of SU{2). First we construct a function g(c) > with g(c) < Vc > c . Recall the 
solution (2.11) 



o(k, t) = <B* At <j cos(At) + i^p-vc - A) sin(At) 



: That is, we just successively concatenate copies of the denning interval along the real axis in order to obtain a periodic function. 
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where A = + c(c — v\)\ ^ 2 and A, v e R, /i e (D, ^ 2 + |/i| 2 = 1. Now denote the explicit dependence on the 
parameters fi, v, A, t as subscripts a^A.t (because a is independent of [i) and compute 



2 

#A,t(c) = ao,o,t(0) - ei xt a 0A A c ) 



sin 2 Wc 2 + + cos(ci) - cos J c 2 + 



(5.7) 



as well as 



A t (c) = Oi, ,t(c) = e 4C *. (5.8) 
We verify for t > and A > that B\ tt (c) = 0(c~ 2 ) for c — > oo, by inspecting both summands separately: 
First, 4c i +X i s hr 2 ^^/c 2 + ^T^j = C( c ~ 2 )i because the modulus of the sine is bounded by 1 and the quotient 

is 0(c~ 2 ) as c — > oo. Second, ^cos(ci) — cos ^^/c 2 + T^^) = ^( c ~ 2 ) because expanding the square root 

inside the argument of the second cosine around 1 cancels to first order the first cosine and the remainder is 
of C(c _1 ) which is squared to yield 0(c~ 2 ). Moreover, for A = 2 and t = 1 we verify £? 2 ,i(c) > as follows: 
Both summands are squares and thus positive semidefinite. The first summand vanishes for c = i/s7r — for 
integer k, while inserting these values for c into the argument of the sine is incompatible with ^ir for integer n. 
Hence £>2,i(c) can be used for the construction in the proof of lemma 3. 

To construct cylindrical functions in LQG whose restriction to homogeneous isotropic connections reduces to 
the functions (5.8) and (5.7) of the homogeneous isotropic connection parameter c, we realize o as the holonomy 
matrix element hu and a as the holonomy matrix element /i 22 in the fundamental representation. Moreover, 
we fix three general points £, v, (, which we denote as superscripts e.g. h^ b , and specify spiral parameters v, A, t, 
which we denote as arguments h^ ab (v, A, t) and define the embedding map r as the linear extension of 

r(l) = 1 (5.9) 

r( e ic *) = /4(l,0,t) (5.10) 

r(B 2A (c)e ict ) = h v n (l,0, t) (h^O, 0, 1) - e^h^O, 2, 1)) (4,(0, 0, 1) - e^fO, 2, 1)) , (5.11) 

where t 7^ in (5.10) and the 1 in (5.9) is understood as the unit function. 

An Explicit Measure for the Symmetry Reduced Sector. Let us now consider the image of r within 
the Loop Quantum Gravity inner product: Observing that r{e lct ) is (upon normalization) a gauge- variant spin 
network function on a different graph for different t, we conclude 

1 



\r(e ict )\ L QG\r(e ics )\ LQG 



{r{e^),r{e^)) LQG =6^-, (5.12) 



where we used \T\lqg = \J (T, T)lqg and the Kronecker delta 8fj\ Moreover, observing that r(B 2 ,i(c) e ltc ) is 
a linear combination of gauge variant spin network functions on a graph always different from r(e ISC ) and for 
s 7^ t on a different graph than r(£?2,i(c) e ISC )), we find 

V(i? 2 , 1 ( C )e itc ),r(e isc )) LQG = (5.13) 



|r(B 2l i(<0 e itc )\LQG\r{e isc )\ L QG 
1 



We have thus established that orthogonality of our generating set and can return to the question of completion: 
For this purpose let us introduce the induced norm / in the span of the generating set: 

H/IU := Hf)\LQG, (5.15) 

which we can freely use to construct a completion of the span of the generating set. However, we have to 
caution that we where not able to prove that the completions in ||.||oo is contained in the completion using 
H-Hind- Hence we can not exclude the possibility that the embedding map r needs to be amended. 



17 



5.3 Asymptotic Analysis and the BKL Picture 

Let us briefly discuss some physical consequences of the construction described in the previous two subsections: 
The BKL-conjecture [5, 6, 7] suggests that the generic evolution of general relativity near a singularity is 
such that spatial derivatives become negligible compared to time derivatives, i.e. each point in space evolves 
separately as a homogeneous cosmology, so the understanding of cosmological models seems to hint towards an 
understanding of the evolution near a generic singularity. Since large time derivatives of the metric imply large 
values of the Ashtekar connection in a fixed trivialization and a fixed chart, a first approach for the investigation 
of the behavior of loop quantum gravity observables near a generic singularity is to investigate the asymptotic 
dependence of observables and in particular spin network functions on homogeneous connections. Thus our 
result, stating that the dependence of a spin network function T(A) on the homogeneous scaling parameter c of 
a connection A = cA* becomes asymptotically almost periodic as c — > oo, suggests that the generic behavior of 
spin network functions near a singularity is well approximated by an almost periodic dependence on the scaling 
parameter c. This scenario has been investigated in the context of LQC [19, 20, 21, 22]. 

However in light of the nonembcddability of standard LQC into LQG one expects modifications to the 
standard LQC dynamics, since it is expected that quantum corrections will become important not at infinite 
momentum scale, but rather at momentum scale of order 1 in Planck units, one expects from dimensional 
analysis 

r(cA,) = V(c)(l + 0(l/c)), (5.16) 

where ip is almost periodic, that the 0(1/ c) corrections are not necessarily negligible when approaching the 
Planck scale and become negligible only when going far beyond Planck scale. It is thus possible that the use of 
definition 3 could lead to observable deviations from results derived using definition 2. However the existence 
of quantum connections of the type defined in equation 5.3 suggests to use definition 3. 

6 Conclusions and Outlook 

Loop quantum cosmology (LQC) is used as a tool for the study of conceptual issues in full loop quantum gravity 
(LQG), because it is a simple toy model that exhibits many features of the full theory. Moreover, since on the 
one hand it is expected that quantum effects of gravity are most significant in the vicinity of general relativity 
(GR) singularities and on the other hand since the BKL conjecture [5, 6, 7] states that the dynamics of GR in the 
vicinity of a generic singularity is well approximated by decoupled homogeneous cosmologies, one expects that 
LQC provides insight into the behavior of LQG where its effects are expected to be most important. However, 
to be able to draw such conclusions, one needs to understand the relation between full LQG and the symmetry 
reduced model. The nature of this relation is not as obvious as one might think, because LQC is obtained as 
a "loop quantization" of a classical symmetry reduction of GR and not as a symmetry reduction of LQG. More 
disturbingly, one finds that the structures of the two theories do not match up: In particular one can show that 
the embedding of the configuration space of LQC into the configuration space of LQG is not continuous [8]. 
Since LQG is a rather restrictive field theoretic framework that does not admit arbitrary changes, we are lead 
to the investigation of symmetry reduced models of LQG in this paper. 

It turns out that the problem of non-cmbcddability of standard LQC into full LQG stems from the re- 
striction to piecewise linear curves in the "loop quantization" of homogeneous cosmology, which in turn implies 
that the holonomy matrix elements, which are the building blocks of the configuration operators in the "loop 
quantization" program are linear combinations of exponential functions of the symmetric connection parameter. 
This motivated our first Ansatz to start with standard LQC, but to allow for small perturbations of the linear 
curves used in the construction. To do this we formally solved the holonomy differential equation as a power 
series in the perturbation parameter, such that the first term in the series is precisely the unperturbed linear 
curve. This provides an intuitive understanding of the relation of standard LQC and its embeddable extension. 
Moreover, it turns out that the formal solution can not converge for large values of the symmetric connection 
parameter unless the perturbation parameter scales with the inverse of the symmetric connection parameter. 

This inverse scaling motivated our second Ansatz, the asymptotic expansion of the solution to the holonomy 
differential equation as an inverse power series in the symmetric connection parameter. The resulting Liouville 
Green series is in general divergent, but one can show that the truncated Liouville Green series captures 
the asymptotic behavior correctly. It follows that the solutions to the holonomy differential equation can be 
characterized as a sum of a continuous periodic function and a continuous function that vanishes at infinity, 
implying that the algebra generated by polynomials thereof can be characterized as continuous asymptotically 
almost periodic function. It turns out that this result generalizes from symmetric connections to arbitrary 
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one-parameter families of connections of form A = cA*^ where spin network functions can be shown to be 
asymptotically almost periodic functions of the parameter c. 

In the following we provided a definition of a symmetric quantum connection that does not tacitly assume 
classical smoothness properties of a quantum connection. The asymptotic almost periodic dependence of confi- 
guration operators on the symmetric connection parameter shows that one has to include at least one nonlinear 
curve in the construction of an embeddable symmetry reduced model of LQG and also allowed us to give an 
explicit construction of an embeddable symmetry reduced model of full LQG, that contains standard LQC and 
a "minimal" extension. This extension consists, in accordance with [1], of configuration observables that vanish 
both, for large and small values of the symmetric connection parameter. 

It is likely that the extension does not change the generic physical behavior in the limit of infinite symmetric 
connection parameter that was established in LQC, because it consists of functions that vanish for large con- 
nection parameters. Thus, elements of the extension are not expected to affect the behavior of the model in the 
limit of infinite connection parameter. Specific physical predictions are however likely to change, because quan- 
tum gravity effects are expected to play an important role not at infinite connection parameter but at values of 
order one (or even less) in Planck units, where the configuration operators introduced through the extension do 
not vanish. In particular, if we are using the usual configuration representation, there are wave functions in the 
extension that vanish for large connection parameter, or physically speaking near the cosmological singularity, 
which are to be expected to yield very different physical behavior than almost periodic wave functions. 

6.1 Outlook 

The results presented in this paper shed new light on the conceptual question what symmetry means at the 
level of the kinematical LQG-Hilbert space Ho and how symmetric states can be constructed thereon. The 
novel answer we get from our analysis is that every specifically chosen symmetric connection cA* induces a 
decomposition of Ho into equivalence classes of cylindrical functions, where two cylindrical functions are called 
equivalent, if their restriction to cA* coincides. Hence we can take the opposite point of view and define a 
particular symmetry by an according decomposition of Ho into equivalence classes of cylindrical functions. 
Starting from this observation there are several future directions for investigation. 

An important task is to examine the mathematical properties of the symmetry-equivalence classes in more 
detail. The ultimate goal here is to work out the general construction of a symmetric Hilbert space "Hsym by 
taking the quotient of Ho and the equivalence relation 12 . It will be important to generalize the construction of 
section 5.2 and to see how a measure on "Hsym can be constructed and to understand, whether the symmetry 
equivalence classes preserve the inductive structure of Ho- 

Secondly, given the explicit embedding of the extended LQC-configuration space into Ho as described in 
section 5.2 we have already proved that the AAP(]R)-functions form an algebra. However one needs to work 
out the action of flux operators of full LQG on a particular embedding of AAP(H) in A. Does this action 
preserve the symmetric sector and moreover the symmetry equivalence classes in "Ho? If it does, the action 
of geometric operators and ultimately the constraint operators on the symmetry-equivalence classes has to be 
analyzed in detail. These questions can be investigated very naturally using the results of this paper, because 
they suggest to define symmetric operators as operators on the LQG Hilbert space as follows: Recall that wc 
have an equivalence relation ~ of cylindrical functions through the pull-back under i : c n- cA*, such that 
any cylindrical function can be written as the sum of a fixed representative in the equivalence class plus a 
function f° whose pull-back under i vanishes. We then call an operator O symmetric if it satisfies 

(Cyh,OCyl 2 ) = ((Cyh + f?),0(Cyl 2 + / 2 °)>, (6.1) 

for all cylindrical functions Cyk and cylindrical functions f° whose pull-back under i vanishes. Notice that this 
relation is by construction linear, meaning that if Oi, 2 are symmetric, then 0\ + a0 2 for a <E C is symmetric 
as is 0\0 2 . 

A third very important goal is to analyze the effect the extension of the configuration space of LQC from 
IR-Bohr to the spectrum A(AAP(1R)) has. Certainly the first goal is to construct a measure on the extended 
Hilbert space. Beside starting from the extended configuration space itself there is the possibility to start from 
our explicit embedding and its usage of the kinematical inner product on Ho using the Ashtekar-Lewandowski 
measure. As described above one might be able to define a measure on H synl and to use it for the extended 
version of LQC. Then one is able to work out the consequences for the physical predictions of LQC. 

12 This generalizes to very general symmetries (e.g. discrete symmetries) and the questions raised here are applicable in these 
cases as well. 
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Lastly, there are some purely mathematical questions: From the viewpoint of differential geometry it will 
be interesting to see, how our procedure (described in section 4.3) for computing the holonomy of a connection 
in the asymptotic regime of a large re-scaling can be generalized to connections on arbitrary bundles. Moreover 
it is necessary to work out its relation to the curvature of the connection as discussed in section C.3 in more 
detail. From a complex analysis viewpoint one is lead to ask when solutions to the holonomy equations are 
analytic in c, this could be investigated e.g. through extending the constructions of section 5.2 to complex c and 
applying results of [17]. From a functional analytic viewpoint one needs to address the convergence properties 
of the embedding map, e.g. through proving bounds that imply convergence. 
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A Boundedness of Formal Solution 

By construction the series (3.30) provides a formal solution to the ODE (3.4). However it is not clear whether 
this series converges, we thus truncate it and investigate the rest term 

e<£ (K,t) := d^^-e^+UM^V^P (A.l) 

k=0 

separately. Insertion into the holonomy ODE yields 



= {J + Ka + fleP + a£>. 

where A^(k, t) = af Kt X^ (t) 2 ^"_i^ . To investigate the asymptotic behavior of e„\ we construct a funda- 
mental system with Wronskian 2 that solves equation (A. 2) to all positive powers in k 

Z (tr) (s,t) = (A.3) 

ia(t) 

with a = ±1, which satisfies 

%) = (K 2 + K a + 7 )Z ((T) , (A.4) 
where 7 (/ S) = leoSpw? ~ 4(«4°ffl) ' Variation of constants of 



yields 



i$ z (+) + ffl = o 



(A.6) 



yielding for c^(k, i) 

e£> («, t) = e£> (k, 0) + J* K(t, r, k) ( x (k, r)e<f> (k, t) + A £> (r)) dr , (A.7) 

where K(t,r,Kj = ^ (Z/ + \(K,t)Zt_\(K,T) — Z/_)(K,t)Z/ + \(K,T)) and x(k,t) = /3(r) — 7(k,t). Imposing the 
initial condition ei £ '(«;,0) = 0, 6u\k, 0) = on the integral equation, we can construct a series solution given 
by — J2m=a ( e n,m+i — e n}mj , where the summands are determined by 

e2>(s,t) = 

- f drK(t,r,K)A^( K ,r) (A.8) 

e n,Li-i(«>*) - en,m(K,<) = J dr K (t, r, n) X {k, t) (eW m (n, r) - e^^K, r)) . 
To find a bound on this series, we impose 

|tf(t,T,/s)| < p(n,t)q{T) 

K := sup < T < t {|g(r)|} 
/C„(/c) := sup < T < t {|p(K,r)g(r)|} (A.9) 

^>(M) := f \&W{K,T)\dT 

Jo 

ip(K,t) := / |x(K,r)|dr, 
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where p, q > 0, such that we can show 

l e L,m+l( K iO — £n,m(K) t)\ 
p(n,t) 

by complete induction. Induction beginning: 



<K<(>W{K,t) 



^iMJ-tfiM)! = |4>)l< / rfr|^(t,T, K )||A(f)( K ,T) 



p(/t,t) / drg(T)|AW(«,T)| 
t 

p(«,t)/C / dr|AW(«,r)| 
Jo 



< P (/t,t)x:^(K,t). 



Induction step: 



(/t,t) - e^,m («,*)! 



p(n,t) 



p{K,t) 



< 



< 



< 



-i— / dr|if(t,r,K)||x(«;,r)||eW m (/c,r) - e^.^^r)) 
K ) Jo 



dTg(T)|x(K,r)||eW m (/t,r) - c^L-i^t)! 
dr q(r) I X («, r) I K ^ («, r) ^° / .r' K^ t) 



(m- 1)! 

< KKZ{k) f dr| X («,r)|^)(«,T)^^ 
Jo (m-lj! 

ll) m - 1 {K,T) 

(m-iy. 



< KKZ{K)<${K,t) [ dr| X («,' 
Jo 

= /c/C(*)^(«,*)^!r^- 



m! 



(A.10) 



(A.ll) 



(A.12) 



where we used the fact that by construction <^(k, and ip(K,t) are non decreasing in the interval [0,t]. The 
series is then bounded by 



ki £) (M) 



m=0 

oo 

< p( K ,t)lC^(K,t)J2 

m=0 



( e ",m+l( K '*) e n,L( K i') y 



ml 



(A.13) 
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Using X(n,t) := info< T <t{|K + ia(r)|}, which is O(k) as |k| — > oo, we can bound K(t,T,n) by 



\K(t,T,K)\ = - (Z (+) («,t)Z ( _ ) (K,T)-Z ( _ ) (/t,t)Z (+) (K,T)) 



K(t-T) + ^ /*«(«)<*« 



-K(t-r)-^/ T *a(a)d« 



< 



v /( K +ia(i))(K+ia(r)) 



(M)V e 
1 

2A(/c,t) 



2X(K,t) 
,(n,t) V 



+ <B~ 



i /, i rt 

< I (P 2 J T 

" 2A(M)v 



|a(s)|ds I 



i/ T *l«(*)l<i* 



I) 
I) 



X(K,i) 

<B2 /o |a(a)|d« j 
= 3-3/0 l«( s )l ds 

A(/c,t) 

Here we have used that k = ic (c e R) is purely imaginary. We can thus choose 

_ cxp(|J^ |a(s)|rfs) 

H^J ~ X(m) 

9 (r) = exp(-i/ r |a(a)|ds) , 

which allows us to set 



K = su Po < r < t {| g ( S )|} = sup < r < t {|(E-^/o T l^)^ 

JC„(«) = su Po < T < f {|p(«,r) g (r)l} = ^[(S'^ ■ 
We thus established the bound 

exp (/ * \a(s)\ds] 



} < exp (±/ * |a(s)|ds 



^(M)|< 



X(k, t) 



Jo |A " )(,t '' 



that ensures that e^(/c,i) — > as |k| — > 00 if a(i), (k, i), %(«;, t) are bounded. 
A.l Estimate of the Error Term 

Let us assume that ip, c are analytic in a domain D(d) that contains all complex numbers z for which min xe /{|^— 
x\ < d}, where / is the integration path for the Liouville- Green expansion and furthermore assume that there 
are constants ki,k 2 ,A s.t. for all x € D(d) we have the bounds 



t) \dr I exp exp 



\f Q Hs)\d t 



Jo\x(s)\ds S 

X(K,t) 



(A.14) 



(A.15) 



(A.16) 



(A.17) 



\i>(x)\<h, [ ds\ip{x)\<k 2 \a{x)\ 

Jo 



(A.18) 



so we can define k = max{fc 1; !f, The bounds yield |V> (m) 0)l < fc if^ and |a> l )(a;)| < through 
Cauchy's integral formula. 

Lemma 4 Given the above bounds, the Liouville- Green coefficients satisfy for t within the integration path and 
n > 2: \air\t)\ < C n ^^±^- with C n = (1 + 2fcd 2 ) Mgjf^ . 

proof: by complete induction following [17]. preparation: 



|di(*)l = ^|/ < \ k 2 



(A.19) 
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Induction beginning: 



\d 2 (t)\ < ljyu\ddu)\mu)\ + \\m\ + liy^Hu)\\d^\u^ 



< 1 (kl + h + Ak 2 ) < | (1 + 2kd 2 ) . 

Furthermore by denoting the digamma function by ip and using Y^k=i k(m-k) = 2 ^+^( m ) as we u 
2 7+V^m) ^ gjj natural numbers m we find 



(A.20) 



as 



i4 m) wi < 



i 



Induction step: 



Moreover, using 



as well as 



we find 



d" 

IF' 



-{d^f 



+ 



l 



^ m \t)\ + 



1| d 



m—1 



-(<*(*)<*! (*)) 



4 <f 

1 m— j. 



2\dt m - 1 

1 I 1 / ™— i N 

< ifc^i + i f 2 |4 m )(t)iidiwi + E (") i4 ro - fe) wiMi fe) wi 

"' 5 x m-l\ (m-/-l)! fci 11 

z=o v 7 

^ jks—(l + k-(l + 2( 7 + i > (m))) + kd 2 
4 a' u \ m 



K+i(t)| < I du\^{u)\\d n {u)\ + \d^\t)\ + \ [ dua{u)d^{u)\ 



2 

< c " (n ~ 2)l i fc^ ^- 1 )' i .u- c " (n - 2)! 

S 2^ 2 2" d»-2 +fe 2« d"" 1 + 2™ d»-2 

< fc ^(!L=i)lf 1 + 2fc rf2 



>n— 1 rln — l 



n-l 



m—1 



2 1 dt" 



TO — 1 



n: 



(to + n — 1)! 



< 



2n+l n 

k (ra + n+1)! C„ 
2 n+1 d m+ ™- 3 n-l 



' /m ~^K*)*»(*)l + I^Wl) < At {m ll~J )l C n (l + k 1,2 



dt 



2^+1 ^jrn+n—l 



n-l 



n-l 



(A.21) 



(A.22) 



(A.23) 



(A.24) 
(A.25) 



If we assume these analyticity conditions and that the initial conditions are such that all integration constants 
Cr'^ in the Liouville-Green recursion relations (3.31) vanish, we can use equation (A. 17), to provide a bound 
on the truncation error that is independent of the Liouville-Green coefficient. This is rather useful, because the 
explicit form, particularly of higher order coefficients, is very complicated. We then get 



< J_ e j>l«MI C B fc£f^. 
~ n n 2 n d n ~ 1 



(A.26) 



Insertion into the truncation error yields the bounded 



e X>H s )l 
£ <0( M )|<_^e*p 



J>lx00l 

\(K,t) 



■ cxp 



*»I«W|]^*^tI*I- 



(A.27) 
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A. 2 Real Analytic Edges 

We assume that components of the connection and of the edge are real analytic in the chosen trivialization. For 
A 7^ 0, we perform a position-independent gauge transformation n i->- g~ 1 T i g , such that m(0) ^ 0. Then there 
are t Q1 e > s.t. inf > T > to {|™(s)|} > e. Specifically, we assume that there exists t x > such that for < r < ti 



p R (r) = Re{\M 2 {r) - ±M(r)) 
/3j(r) = 7m(±M 2 (r) - |M(r)) 
ai? (r) = Re(n(r) - M{t)u(t)) 
oli{t) = 7m(n(r) — M(r)n(T)) 
are real analytic for > r > t\. Then for each ieN there exist 

> o s.t. > suPo < T < t , H {|^ii)| } 

4,4, > s.t. a) >su Po < T < t:vj {|^^|} 
fitJpR > s.t. fa R >sup < r < t , R {|%^|} 
0},%, > s.t. fa\ > SU p < T < t , 7 {|^ll|}, 
which lets us define on := max{a^, a}}, fa :— ma,x{fa R , faj} and 

,o ,o 

We thus have the bounds 



tk ■— min{i° K , t° a j, t°p Rl t° jI} t„ R , t k aI , tp Rl tpj, t , t\}. 



which implies 



13 



\a(t)\ < a for t < t 
\fat)\ < fa for t < t , 

\ X (K,t)\ = \fat)-l(^t)\<\P(t)\ + \ 7 (K,t)\ 

a 2 (t) 3d 2 (t) a(t) 



< fa + 



16(k 



2 J 



4( K +^) 



a2 



1 2 3 aj 
- Po + -^o + j6A 2 (K,t) + 4A(/s, i) 



for t <t a 
for t < i 2 . 



|AW(«,t)| - 2|<e!/o«(^if(t) 

(3 i 0) 2 |e*/o* - ±a 2 (i) - fd(t)) 

< 2( B ^ t (^ + ia 2 + ia 1 ) fort<ti. 

= / \A[ e) ( K ,T)\dT 
Jo 

< J" + \o? + \a^J (ce^ 4 - l) for t < h. 



13 Again, \(k, t) := inf < T < t {|/{ + ^a(r)\}. 



(A.28) 



(A.29) 



(A.30) 



(A.31) 



(A.32) 



(A.33) 



(A.34) 
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We thus have a bound on the truncation error 
l £l l K '^l^ e V s U 2a„AfK.fl 6XP 



€3«ot . (iQ{3 o \2( Ki t) + 4a 2 \ 2 ( K ,t) + 3aj + 4a 2 \(K,t)) ^ 



A 3 ( K ,i) 

(A.35) 

valid for t < t 2 which, using that by construction A = C(|c|)| for c — > oo = |«| — > oo, establishes 

\ef } {K,t)\ -s-0 for |c| -> oo. (A.36) 

B Perturbation Approach 

Restricting to a particular type of edge refers to a background. Even worse it prevents LQC to be continuously 
embedded into LQG. For this, one would need the solution to (2.7) for the isotropic case but for arbitrary edges. 

In this section we are now going to construct the general solution to (2.7) in terms of a perturbation 
around a known explicit solution of section 2.2. For this we use methods of [14], similar to the Liouville-Green 
approximation described there. However instead of giving an approximate solution we will compute the exact 
solution in terms of a formal power series in the perturbation. 

B.l General Perturbation 

Suppose a solution to (2.7) is given for a special edge eo(t) — (eo 1 (t), eo 2 (i), eo 3 (t)) =: (xo,yo,zo). Now 
assume the edge e is deformed into another edge e such that e(t) = (x + ex, y + sy, z + ez), where 
x = x(t), y = y(t), z — z(t), and e = const is a deformation parameter. Then we can write the components of 
the tangent vector e as 

x = x + e x y = yo + e y z = z + e z . 

This implies 



m = (x - iy ) + e (x - iy) n = z + e z 

=: m + e m =: n n + e n 

Moreover we introduce 

M := ^ and M := ^ with M = A _ MM 

u mo mo mo u 

M = A-2M (M + MM ) 



(B.l) 



in order to expand 



m rh + em 1 , ^, 1 

M = — = — = (mo + em) 



m m + em m Q y \ _ (_ e M) 

oo 

(M + e(M + MM )) ^(-l) fe e fe M fe 



fe=0 



CO 

= ^(-l) fe {e' £ MoM fc + e k+1 (MM k + M M fe+1 )} 

k=0 

oo 

= M + ^(-l) k e k ^M)M* r -MM k - 1 -MtfM*} 



fe=i 

M °° ~ 
= M - ^V(-l) fe £ fe M fe 

where we need to assume \eM\ < 1 from the first to second line in order to ensure that the geometric series 
converges. Notice that this assumption is re-parametrization independent: If we introduce a new edge parameter 



2G 



T = T(t), then for an arbitrary function f(t) it holds that f(t) = ^f- = d ^jp T(t), hence the T-terms cancel 
by construction of M. As a consequence, it is straight forward to write down 



M 2 = M 2 -2M "^(-l)VM k + -t ^^(-l) fc+; e fc +<M fc+i 
M k=i LMJ k=ll=1 

M °° ~ riwi 2 00 ~ 

= M 2 -2M ^ s T(-l) k e k M k + ^ V(r- l)(-l)V"M r 

MtC LmJ ^ v 

= M 2 + (2M M) £ + ^ V(-l) r e r M r ((r - 1)^ - 2M 
M^ I M 



and 



~~ - i - 2 00 ~ 00 

M = Af - MM Z M Y{-l) k e k M k - ^ V (-l) fe e fe fc M^M 



M 2 



fe=i 



/ M 


"M" 


") 


VM 


.M. 





fc=l 



^(— l) fe e fe M fc - |£ ^(-l) fc e fc fc M fc 



2 00 



fe=i 



^2 



= M + M e-M" 2 ^](-l) fe e fc M fe {(fc- 1)M + MM} 

fc=2 

Finally we use the expansion 

(x + e x) 2 + (y +ey) 2 + (i +ez) 2 



\m\ 2 + n 2 



2/0+^0 + 2(x a; + yay + z z) e + (x + y + f") e 2 



=: \m Q \ 2 +nl + 



2Y 



e + (\m\ 2 +n 2 ) e 2 



in order to express 



N = 



ic(n — Mn + ic(mm + n 2 )) 
= -c 2 (\m \ 2 +nl + 2Ys + (|m| 2 + n 2 )e 2 ) + ic(n + e n) - ic(n + e n) M 

= -c 2 ^. • • ) + ic(n +s n) - ic(n a + e n) |m - ^ ^(-l) fe e fe M fc j 



fc=i 



= -c 2 (\..) + ic(n + e n-M (n + en)) + ic^£ ^(-l) fc M fc {n e fc + n e k+1 ) 

M k=i 



= -c 2 (\..) +ic(n + e n-M (n + en)) + ic( - n M e + ^= (n - ^(-l) fc M fc e fc ) 



M 



fc=2 



-c 2 (|m | 2 + Uq) + ic(n - M n ) + { - 2c 2 F + ic(n - M n - n M)} e 

M 



A/7" ~ 

+ { - c 2 (\m\ 2 + to 2 ) + ic M(n M - n) } e 2 + ic ~ (n - ^(-l) fc M fe e 



fe=3 
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So we are finally able to expand the bracket term in (2.7) into orders of e : 
jlM 2 -iM + iv} = ^M 2 -^M - C 2 (\m \ 2 +n 2 )+ic(h -M a n ) 

+£ { \ (M M — M) — 2c 2 Y + ic(n - M n - n M) } 

+e 2 |i((M) 2 - 2M Mm) +i((M) 2 + MA?) - c 2 (|m| 2 + n 2 ) + ic M(n M-H)} 
+ ^(-l) fc M fc e fe j-MM^^fc-ljMM- 1 - 2M ) + -M~ 2 ^(fc — 1)(M) 2 + MM^j 

k=3 

+ icMM- 1 (n -nM- 1 ^ , (B.2) 
where we have assumed \eM\ < 1 and Y := xqx + yoy + z z. 
B.2 Simplifications 

As it stands (B.2) holds in full generality. However, in order to solve (2.7) with the perturbation Ansatz, we 
require Mq and no to be constants. Hence Mq = = no- Moreover we may assume that we perturb eo{t) only 
in perpendicular direction, that is for the edge tangents eo,e it holds that (eo,e) = (eo,eo) = |eo | 2 - Then it 
follows that Y := x^x + y^y + z z = 0. 

If we leave the beginning point of the edge fixed under perturbation we have e (t = 0) = e(t = 0). More- 
over without loss of generality we can choose e (t) such that e (t = 0) = e(i = 0). 

Finally we are free to choose a convenient parametrization, for example arc-length parametrization of e (t), 
such that | mo | 2 + n§ = 1. We will refrain from the latter and write using 



k := ic 



{^M 2 -^M + N} = ^M 2 + K 2 (\m \ 2 +n 2 ) -nM n 

= +e |-(M M - M) + n(n - M n - n M)| 

+e 2 |i((M) 2 -2M MM) +^(M) 2 + Mm) + n 2 {\m\ 2 + n 2 ) + k M(n M - n) } 

+ ^(-l) fe M fc e k {-MM" 1 ((fc-l)MM^ 1 -2Mo) + -M" 2 ((fc - 1)(M) 2 + Mm) 

fe=3 

+ KMM-^no-nM- 1 ^ 

OO 

=: /C 2 ( K ) + 5> fc / fe (M) , (B.3) 
fc=i 

where we have introduced the obvious shorthand in the last line. In particular 

IC 2 (k) := ^M 2 + K 2 (\m \ 2 + n 2 )-KMo n (B.4) 

B.3 Solution by the Liouville-Green Method 

With the expansion (B.3) we can write (2.7) as 

oo 

d= {/C 2 ( K ) + ]T£ fe /*(«,*)} d (B.5) 
fe=i 



28 



Now we make a Liouville- Green Ansatz, similar to section 3. That is we set d — d(n,t) = E d^\n,t) with 

<T=± 

CO 

a— ± a— ± n— 

oo 

(T— ± tr— ± n— 

oo 

d{n, t)=J2 ^ («. i) := X] eff0 E d ™ CT) (M) + 2 tr £ («. f ) + d { n ] (k, i)) • e" (B.6) 



cr=± (T— ± n— 



In the last line we have used the fact that a 2 = 1. As we will see at the end of this computation the 
choice of initial conditions 

d( CT )( K; 0) = l and rf (CT) (K,0) = , (B.7) 

which have to hold for arbitrary e, will be convenient. Now we plug the Ansatz (B.6) into (B.5). For clarity we 
will again suppress the dependence on t, k and simply write e.g. d^ instead of d^\n, t). This gives: 

oo oo 

E ^ £ {K^r+2aKd^ + #) -J&ifr- d^^^Y = 

cr=± n=0 fe=l 

This has to hold at any point of e(t), that is for arbitrary values of t. Hence, it must hold that 

oo oo oo n 

£ {2 * * + ^ E a ^}^ = E £ "H a)+2<T/c - E } + + 2 * jc 4 ff) = o 

n=0 fe=l n=l fe=l 

This has to hold for arbitrary values of e, hence separately in ever?/ order n of e: 

Q(e°)\ + 2 <T /C = (B.8) 

n 

fc=i 



g(g" >0 ) + 2 ^ AC = A<f> with := E €l k A (B-9) 



B.3.1 Solution to O(e ): 

This is a homogeneous linear ODE of second order with constant coefficients. Its solution is given by 

with integration constants Aq <t ' ) , . 
B.3.2 Solution to O(e n>0 ): 

This second order ODE is inhomogeneous but still linear with constant coefficients. Its general solution d^ 
can be obtained as a linear combination 

4 ct) =<^om + < S p (B.10) 

of the general solution c^jom to tnc homogeneous equation plus a special solution d^ SP to the inhomogeneous 
equation. The homogeneous part is equivalent to the O(e )-case and given by 

A*) _ a(<?) , dW (p -2<7/Ct 

n,HOM — n.HOM ^ n,HOM * 

with integration constants ^^"hom' ^nHOM' The special solution d^ SP can be obtained from d^ HQM by the 
method of variation of constants, that is we make the Ansatz 

di% («, *) = («, t) + («, i) e- 2 ^ 4 . 
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With the usual requirement: \T\ A ( n ] (k, t) + B ( n } (k, t) e = this gives 

<*£spM) = (-2 < 7/CBi ff) ( K ,t)+4rf( K ,())e- M 

If we plug this Ansatz into (B.9) we get (recall that a 2 = 1): IT_ B$?\n,t) = -jhc^^ ^n\n,t) and thus 
in turn 



where A^$ p , B^ p are again integration constants. Therefore the special solution d^ SP (n,t) is given by 

41p («> = <SP + B Sp + ^ / (1 - ^ * («, -) d S 

Using (B.10) the final solution of (B.9) is given by 

4r' = + sSom + ^ ( (1 - AW ( k , ,) ds (B.ll) 

where we have introduced the overall integration constant D^g P := A^l iOM + A^g P + B^g P . 

B.3.3 Implementation of Initial Conditions 

Now we are going to implement the initial conditions (B.7) into (B.ll). 

oo oo 



d<"> («, 0) = 1 = VJ W («. 0) • e" = A a) + + E «sp + ^Aom) • ^ 

n— 71=1 

OO OO 

d<->(«,0) = = ^dW( K ,0).e»,-2^(5W + ^BW OM . e ») 



n— n=l 

As these conditions have to hold for arbitrary e, it follows that 



A) - 1 ^0 — U anCl n.HOM - U ~ -^n.SP 



Therefore 



<^)( /S) i)= e *Ktjl + Jl^ e » J* ( 1 _ e 2*K(.-t)) \W( K , S ) d*| (B.12) 

n 

where A£^ (k, t) := d^]_ k (K, t) t). Moreover by construction («, i) = 1 and is given according 



fc— 1 

to (B.3) and JC is defined in (B.4). 
B.4 Solution 

Now we are set up to construct the solution to (2.4) for an arbitrary edge e(t) using the Ansatz (2.8) in terms 
of a perturbation about an edge eo(t). According to section B.2 we demand 

e (0) = e(0) and e (0) = e(0) . 

This implies m(0) = mo(0). Imposing the initial conditions (2.5) we obtain: 



o(k,0) = 1 = A (+) df- + \n,0) + A ( _)d(-)(K,0) = V™o(0) (A {+) + A { _ } ) 
a(0) = Kn(0) = A (+) d (+) (fi;,0) + A ( _ ) (i ( " ) (K,0) =K,{A {+) - 
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Hence Ai a \ = i I" , 1 + q- ^" ^) ] anc j we g e t the final solution 

1 + ^E e "7 (l-^^A^s)^ (B.13) 

71=1 ' 

with \^n\n, s) and /C = JC(k) given below (B.3) and m given in (B.l). 
B.5 Perturbation about a Line 

Certainly the functions fk(n,t) of (B.3), which are needed in order to explicitly compute (B.13), are still quite 
complicated. Also the definition (B.4) of JC involves a square root taken from a complex number. To simplify 
this situation we can construct a solution of (2.7) for a general edge e(t) as follows. Given e(t) we construct the 
solution to (2.7) as a perturbation of a line eo(t), for which again 

e (0) = e(0) and e (0) = e(0) (B.14) 

holds. W.l.o.g. we can choose the maximal simplification for e being a line. That is 

eo(t) = (x {t),0,0) e(t) = (x (t) , ey(t) , «(<)) (B.15) 

Moreover we choose arc-length parametrization of the line, that is we demand 

1 = ±1 = m(0) = m (0) (B.16) 

Then we have 



tuq = xq = 1 Mq = = mo "-o = = no K- 2 = k 2 = —c 2 m = y 
rh = y M = m = y n = z 

Under these assumptions, expression (B.3) can be simplified to 

1 _ _<) 1 , ' , _ "1 o <-) ( 3 , n 1 



3,, _ w ~o It 

e" (-1)" \ 

fe=3 

and expression (B.13) reads 

oo rt 



Notice, that if we set 



(B.17) 



{-M 2 --M + 7vj = K 2 + e n?i + e 2 | -(M) 2 + -MM + n 2 (\rh\ 2 + n 2 ) — n Mn| 

+ Y, £k (-^M^ 2 {-(k - l)(M) 2 + -MM - k Mn} (B.18) 



^^E^^Il + fE 6 " Al-e 2 ™ (-*)) AW(«,s)ds} (B.19) 



then for c — > oo the solution (B.19) takes the form 

o(/s,t) = ^ (E ,TK *{l + C'(c- 1 )} . (B.20) 

cr=± 

In fact, this property holds for arbitrary curves in the limit c — > oo. The proof is given in section 4. 

C Geometric Interpretation of the Parameter c 
C.l The Ashtekar Connection Revisited 

For completeness we will briefly sketch some of the insights obtained in [15], where a coordinate free treatment 
of the Ashtekar connection is developed. Notice that unlike in the rest of this paper we denote vectors by the 
symbol V. 
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We are working with a (3+l)-decomposition of spacetimc (A4,g) = Ex S. The Cauchy surfaces £ are ori- 
entable. Hence at every m G M we can decompose the tangent space T m A4 = T m A4^ © T m M.^ into two 
orthogonal subspaces, the component T m Ai" = T" m S tangent to £ and the normal component T m A^^ = JVn, 
where n is the surface surface normal unit vector and N G 1R. With the signature (s, +, +, +) of Ai (s = ±1) 
we have g(n,n) = s and g(n, X)=0, g(X,X) > for every X g T m £. To shorten notation, we use the 
isomorphism between the local tangent spaces T m S ~ R 3 and the Lie algebra so(3) induced by the equiva- 
lence of the defining representation of 5*0(3) on 1R 3 and the adjoint representation of 5*0(3) on so(3). Let 
{e/}/ = i.2.3 be an orthonormal basis-frame of T m E, e.g. the Cartesian standard basis of R 3 . Moreover let 
be a basis of so(3) <~ su(2), orthonormal with respect to the Cartan-Killing-metric thereon, e.g. 
(ti)jk = —cijk- Now we choose a fixed identification ej -H- 77 and can write any element X € T m S equivalcntly 
as X = ^2 j Xiei o Y^i Xiti, where Xj = g(X, ej) denotes the expansion of X into {ej}j=i,2,3- Consequently 
we will write capital Latin indices to denote both so(3) indices and indices with respect to the orthonormal 
frame {ej}. Then we can decompose the covariant derivative V M on M coming from the Levi-Civita-connection 
as 

V^e, = (V£e J )" + (V£e J )- L 

= E^ e ^^) e J + sg(V%e I ,n)n 
J 

= ^2,9{^ K ei,ej) ej + s g{K{e K ,e I ),n) n , (C.l) 
j 

where in the last line we denote by V s the covariant derivative coming from Levi-Civita-connection on E, which 
is for any ej, ej € T m T, precisely given by the tangential component of V M . Now we introduce the shorthands 

k K i = k{e K ,ei) = g{K{e K , e/),n) := ff(V£e/,n) = -g{e u V^n) =: -g{e u W n (e K )) , (C.2) 

where Uki are the components of extrinsic curvature and W n : T m S — > T m £ is called Weingarten map 14 . Notice 
the symmetry of km = \{kxi + kix) = k^xi)'- By construction torsion T M of the Levi- Civita- Connection W M 
vanishes and we have 

T M (e K , ej) = = V^e, - V^e K - [e K , ei ] 

and therefore kxi — kjK if the vectors e^,e/ commute, that is [ejf,e/] = 0. For the Ashtekar connection we 
have for any two smooth vector fields X, Y e r(T£): 

Vx^ = V|F + /3 W n {X)»Y (C.3) 

where /? is the Barbero-Immirzi parameter and we have introduced the product: 

X . Y = XiYj ej • e j = ]T X i Y J ^JK e K , (C.4) 

7, J 7,J,7f 

where we have expanded X = X^75(^ e 7) e 7 — : ^7^7 e i- We can write this for the orthonormal basis 
{e/}/ = i i2 ,3 as 

7, 

= V^ej - /3 51 fc 7L £ L jm e M (C.5) 

L,M 

or in components using g(ej, ej) = <5/j 

A/jjf := (^7)j7f =g{V^ej,e K ) = T IJK - P^k IL e LJK ■ (C.6) 

L 

Here T^/j correspond to the usual Christoffel symbols 15 . 

14 By construction W n has only components in T m £. We have g(n,n) = s = const and hence (V££ <?)(n, n) = = 
ex(s) +2g(V^n, n), where we use the symmetry of 3 and the fact that V-^ is metric. 

=0 

15 Its components can be obtained for an arbitrary basis from the well known Koszul formula [23]: 
Tijk = g{S%ej,e K ) = §( eig(e.j,e K ) - <^K ff(ej,ej) + e J g(e K ,e I ) -g(e/, [ej,e x ]) + 9(e/r , [e/,ej]) + s(ej, [e K ,e/]) ) 



32 



C.2 Curvature of the Ashtekar Connection 

Our previous considerations imply for the curvature of the Ashtekar connection 
R A {e!,ej)e K = V^V^e* - V£v£e* - V£, iej] e* 

= ^ {(ei(AjKM) — £j(Aikm)) + ^ (AjklAilm — AiklAjlm — CulAlkm)} £m 

M L 

= R{ei,ej)e K + P^{{ei{k JL ) - ej(k IL ))e KL M + ^{Tij L -Tj IL )k LN £ KNM } e M 

L,M N 

+/3 2 ^ (km k JM - kiM k JK ) e M (C.7) 

M 

where Cijl = ^ijl — ^jil = 9([e/, e.j],eK) from the torsion freeness of V s . This gives for the scalar curvature 
of V- 4 

R A := JX^e/.ejJej.ej) 

= H { M A ^r) - e ^(^.//)) + ( A JJlAili - A IJL A JLI - C IJL A LJI ) } 

I, J L 
IJ 

= i?-/3 2 {(tr(fc)) 2 -tr(fc 2 )} , (C.8) 

where the term proportional to (3 vanishes, because it is antisymmetric in /, J. This can be compared to the 
decomposition of the 4-dimensional curvature scalar under a 3 + 1-decomposition M = E x R of the manifold 
[15]: 

R M = R z + s ^ tr ( k )Y -tr(/c 2 )) +2sdiv(-tr(fe)n + V^n) 

= i? A + (s + /3 2 )|(tr(fc)) 2 -tr(/c 2 )| + 2sdiv(-tr(/c)n + Vfn) . (C.9) 

Here div(X) = J2i v(^^ X ' e i) + s di^^X, n ) denotes the divergence of a vector field X e Y{TM). Starting 
with the Hamilton constraint given in [2], p 124 as 

H = --JL=(K l a K> KiK\)E«E\ JteTqR (CIO) 

where I = 1,2,3 denote sw(2)-indices and a,b denote indices in TE, moreover the densitized triads Ef are 
given in terms of the non-densitized inverse triads ef as Ef = Vdct q ef. This can be rewritten as 



H 



Vdct q 

Using (C.8) this can be written as the well known form 

/? 2 



s(tr(^ 2 ) -(ir(if)) 2 ) -R (C.ll) 



Vdct q 



H = sR A + (f3 2 - s)R . (C.12) 



This finishes the short digression on the coordinate free description of the Ashtekar variables. We will use this 
formalism in order to re-obtain the geometric meaning of the parameter c in case of the homogeneous isotropic 
Bianchi I -universe. 



C.3 Interpretation of the Symmetric Connection for Spatially Homogeneous Isotropic 
Models 

In what follows we will examine the properties of the Ashtekar-connection on Friedmann-Robertson-Walker- 
spacetimes, whose metric tensor can be written on Mi = I x E (I c It, open) as 

g = s dt ® dt + a(t) 2 gv , (C.13) 
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where (E,<7s) is a connected Riemannian manifold with constant sectional curvature Ky. — k and gs is inde- 
pendent of t. At every m € A4i we decompose T m .A/f / into the tangent space T m T, and an orthogonal direction 
spanned by the surface normal n = Assuming [n, A] = V A e T m £, we have 

g(W n (X),Y) (C = 2) ff (V^n,r) 

footnote 15 1 9 



--( a( t)^(A,r)) 



Hence (C.3) reads 



a(t) a(t) gx(X,Y) ~» W„ = ^ Id TS . (C.14) 



v£y = v^y + p - x»y (c.i5) 



or in components as in (C.5) 

A IJK ■= (Ai)jk = ff(V^ej, e K ) = Tijk - k IL e LJK 

L 

= \{Ckji + C IJK - Cikj) -P~ a Y. ^lcljk (C.16) 

L 

where we have used the product structure e/ • ej = J2i j e UK &k- Moreover Cijk — g([ei, ej], e/f), and the 
result for T IJK follows directly from footnote 15. In case of Bianchi I we have T IJK — and g(d m ,d n ) = a 2 5 mn 
hence the coordinate basis {d m } can be related to an orthonormal basis {e/} of TT, by d m — J2 m e mi e i = 
aS mI e, . Hence A IJK = -P%Y,l^ljk = ( A fr M )jK = -Afe MJ K- Therefore Af = and 
A b Ej e biAf =/3 a f E/ = /? a <5f. Therefore we get the usual [13] identification 

c = l3a (C.17) 

This yields for (C.7) in an orthonormal frame {ej} on TS 

R A (e I ,e J )e K ^R j: (e I ,ej)e K +l3 2 Q (a • ej) • e K (C.18) 

and we consistently obtain for Bianchi I with i? E = (C.8) 

R A = -6/3 2 ^ & c 2 = i?- 4 (C = 8) ^{(tr(fc)) 2 -tr(fc 2 )} = (/3a) 2 (C.19) 

where we have used the fact that kn, — —^il an d dimE = 3. Evaluating (C.9) for Bianchi I we find in 
agreement to [24], p 97: 

*" —•((:)'+!)• < c - 2 »» 

This shows, that in a situation where a <C a one can understand the limit c — > — oo as a blow up of scalar 
4-curvature on M. , that is a situation close to a singularity. The fact that in the limit c->oowe re-obtain the 
LQC framework indicates an affirmation of the BKL-picture, as discussed in section 5.3 . 
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